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Abstract 

Boundary conditions for unconventional superconducting order parameter 
are obtained on the basis of a microscopic theory. The upper critical field in 
a superconducting film is examined for unconventional superconductors with 
two-component order parameter and is compared with the case of accidental 
degeneracy. It is shown, that for both cases temperature dependences of the 
upper critical field in a superconducting film may depend substantially on 
the quality of boundaries, and under certain conditions have a kink due to 
the influence of film boundaries. The location of a kink point appears to be 
dependent on a film thickness. If for a massive sample there is another reason 
for the existance of a kink, the interplay of the reasons may lead to a specific 
behaviour of a kink point location as a function of a temperature or a film 
thickness. A new test is suggested permitting to distinguish between E\ and 
Ei types of pairing in a hexagonal superconductor near T c . 

PACS numbers: 74.20.De, 74.60.Ec, 74.50.+r, 74.70.Tx 
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I. INTRODUCTION 



The properties of unconventional superconductors are currently widely discussed. In 
particular, the types of pairing which correspond to two-component superconducting order 
parameter rj = (771,772) are of interest especially due to the unconventional properties of 
some heavy-fermion superconducting compounds exemplified by UPt^ . Another possibility 
of treating the properties of these compounds is related with the accidental degeneracy, when 
two different types of pairing (with one-component order parameters 7/1 and 772 correspond- 
ingly) have for some reasons very close critical temperatures T cl , T c2 (see, for example, |]-§J]). 
Splitting of the superconducting phase transition |||6|], a kink (a change in slope) in H c2 (T) 
and H c i(T) []7 HlO|1 , several phase transitions within a mixed state [|7|,p~l1 [TT|1 or under the 
pressure [|TB[ are the important experimental evidences in favour of unconventional pairing 
in UPt^. They can be described to some extent within the framework of several alternative 
approaches |T||2l,|iT)|-p5|] . Though the existence of a tetracritical point on H-T -phase diagram 
of superconducting UPt 3 at any orientation of a magnetic field leads to failure if one bases on 
a theory with two-component order parameter, the possibility of two neighbouring critical 
points (instead of one tetracritical point) is compatible with the theory P,[7|,|8| ri5|ri7|JT8[] . So, 
for the unambiguous identification of the type of pairing in UPt^ additional experiments are 
needed. 

In this paper we examine a temperature dependence of the upper critical field in a 
thin superconducting film, when hexagonal or tetragonal superconductor is supposed to be 
unconventional with two-component order parameter or with accidentally nearly degenerate 
order parameters. It will be shown that the influence of film boundaries themselves may 
result in appearance of a kink in H c2 (T) at a temperature T = and the conditions 
under which it occurs will be determined. The value of depends upon a film thickness. 
For T cl 7^ T c2 the kink actually takes place in a massive homogeneous sample as well. 
The behaviour of the kink in a film in comparison with the one in a homogeneous sample 
has some distinctive features. Massive superconducting samples of UPt^ appear in fact 
to be nonhomogeneous, containing a macroscopic incommensurate structural modulation 
2T| . Considering a thin film as the simplest example for examining the influence of a 



nonhomogeneity upon H c2 (T), one may come to the conclusion about a possibility of a 
noticeable contribution of lattice superstructures in UPts to the observable kink of H c2 {T). 
The presence or absense in experiment of the determined below temperature dependence of 
H c2 for thin superconducting film of UPt 3 (or of any other superconductor) could help to 
identify the type of pairing in the sample. 

For a solution of the problem it is nessesary to make use of the boundary conditions for 
the order parameter in unconventional superconductors. This problem is examined in the 
second section of the paper. For this purpose the microscopic approach is used, based on 
the quasiclassical Eilenberger equations. 
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II. BOUNDARY CONDITIONS FOR THE ORDER PARAMETER IN 
HEXAGONAL UNCONVENTIONAL SUPERCONDUCTORS 



The Ginzburg-Landau free energy for a hexagonal unconventional superconductor with 
strong spin-orbit coupling and with two- component order parameter (7/1,772) can be written 
up to the second order invariants as follows |30|JT|] : 



F = J dV(arjir]* + K x p\r])p^j + K 2 p*r]*p j r ]j + K 3 p*r]*Pjr]i + K 4 p* z r]*p z r]i). (1) 

Here a = a(T — T c ), pi = —ihd Xi — (2e/c)Ai. The z axis is directed along the hexagonal crys- 
talline axis. The indices i,j are equal to 1 and 2 and correspond to the x and y components 
for the operator p. We let H = 1 in this Section. 

In the case of one-component order parameter the corresponding free energy may be 
obtained from (1) by setting r/ 2 = and K 2 = K 3 = 0. 

The Ginzburg-Landau equations stemming from (1) should be supplemented by the 
boundary conditions for the order parameter. These boundary conditions must be consistent 
with the absence of the current normal to the boundary at the boundary with vacuum. Since 
j = —c5F/5A, the expression for the superconducting current has the structure: 

j = -2ie(77*Jii + if 2 h 2 - c.c), (2) 

where hi, h 2 are linear combinations of dtfj. The condition of the vanishing of the current 
normal to the boundary can be written as: 

f hn = bnVi + bi2V2 ,^ 
\ h 2n = hiVi + h2V2 

The quantities in (3) are taken at the boundary, subscript n denotes the projection on the 
normal to the boundary n which is directed inwards for defmiteness, 611,622)^12 = 621 are 
real. The boundary conditions (^) are equivalent to ascribing of the effective surface energy 
to (1) with density 6^77*77^. Within the Ginzburg-Landau approach 6y are independent upon 
temperature and they are functions of relative orientation of the crystalline axes and n. The 
values of 6^ must be derived from microscopic equations. This can be done without loss 
of generality, for instance, in the absence of the external magnetic field and for the order 
parameter varying only in the perpendicular to the plane boundary direction. It permits to 
decouple two equations for the two-component order parameter, which substantially simpli- 
fies the consideration. Let's denote by I the coordinate along perpendicular to the boundary 
direction and take I > for defmiteness. The generalization on the case of the presence of 
magnetic field and of spatial variation of the order parameter parallel to the boundary is 
obvious due to (3). Generally speaking, in the absence of the external magnetic field the 
current parallel to the boundary may flow, so the magnetic field may be induced. Naturally 
near T c the account of the internal magnetic field results in corrections to the boundary 
conditions which are small in the measure of (T c — T)/T c . So such effects usually may be 
neglected and we don't take them into account. Below we shall also consider for simplicity 
the case of the spherical Fermi-surface. 
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When the pairing is singlet, Eilenberger equations for anisotropic superconductors ( see 



e.g. [^1|,0]) can be reduced to the following system: 

' {2uj m + v F V R )f(p, R, uj m ) - 2A(p, R)g(p, R, uj m ) = 
< {2Lu m -v F V R )f + (p,R,u; m )-2A*(p,R)g(p,R,tj m ) = 
, v F V R g(p, R, uj m ) + A(p, R)f + (p, R, uj m ) - A*(p, R)f(p, R,oo r , 



(4) 







Here u m = (2m + l)irT is the Matsubara frequency, p = p/\p\ is the unit vector in the direc- 
tion of Fermi momentum, v F (p) is the Fermi velocity, A(p, R) is the gap function. Anoma- 
lous Green functions f(p,R,u m ) and f + (p, R, uj m ) = /*(— p, R,u m ) and Green function 
<?(p, -R, cj m ) = g*(—p, R, oj m ) satisfy the normalization condition 



g 2 (p,R, uj m ) + f(p, R, u m )f + (p,R,uj r , 



The equation of self-consistency is: 

A(p,R) = -nTY: 



d 2 S' 



(2tt) 3 v f 



(5) 



(6) 



where V(p,p') is the anisotropic pairing potential and the integration is carried out over 
the Fermi surface. The system (^) should be supplemented by boundary conditions for 
quasiclassical propagators p3|-|35|| . If the nonmagnetic impenetrable boundary is specularly 
reflecting, these boundary conditions take the form: 



g(p) = g(P), f(p) = f(p), f + (p) = f + (P)- (7) 

Here p and p denote the incident and reflected momenta correspondingly. According to the 
mentioned above, expressions for bij may be derived without loss of generality in the case 
of the real gap function A(p,R). Then it is natural to introduce fi = (f + / + )/2, fa = 
(/ — f + )/2. Neglecting higher powers of A near T c , we come from (4) to the following 
equation: 

/i(fcJ,w m ) - ^dfh{p^u m ) - = 0, (8) 

where v\ is the /-component of the Fermi velocity v F . 

The solution of this equation satisfying the boundary conditions (7) is: 

A(p,/,^ m ) = — / {ex P (-|^i|Z-Z'|)A(p,0 + exp(-|^|(/ + 0)A(p,OK. (9) 

\Vi\ Jo Vi Vi 

In the case of one-component order parameter the gap function is supposed to have 
a factorized form A(p, R) = i](R)i/j(p), and for the two-component order parameter the 
corresponding form is A(p, R) = r)i(R)ipi(p) + 772 (i^V^p)- Thus the symmetry of the gap 
function in homogeneous superconductor is defined by the symmetry of ipi^ip)- Substituting 
(9) into (6) we get with account of these representations two integral equations in the case 
of two-component order parameter: 
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n\ 7rTA f°° f rfa i /m/ A (P' Z/ ) i i 2 ^ m n/ v\\ , 

*(0 = J^p^S ^L J d SMp){^ r eM-\ — \\l-l\) + 

+ expH^|(Z + /'))}d/'. (10) 

Here A is the effective coupling constant in the anisotropic case: 

mp) = -J v (p>p'^Wj0^- ( n ) 

It is related to the critical temperature T c by irT c \J2 m |^m| _1 = 1- The linearized equations 
(10) are applicable at the distances less than £(T) from the boundary. The function f 2 does 
not contribute to them being odd in u m . 

For the case of the one-component order parameter the equation (10) and corresponding 
boundary conditions for pi — Py-type of paring in tetragonal superconductor have been 



recently examined in fl36|. For those orientations of the boundary for which the relation 
■0(p) = ~'4 , {P) holds, the solution of (10) in the case of the one-component order parameter 
is r]{l) = CI. Similarly if the order parameter has the property ip(p) = tp{p) the solution of 
(10) is T)(l) = C. The momenta of electrons reflected from the boundary are related to the 
incident momenta by the reflection accross the plane parallel to the boundary. Since there are 
seven planes of symmetry in the symmetry group Dgh, we may write the boundary conditions 
using the characters of the reflections accross these planes (for seven corresponding crystal 
orientations with the normal to the boundary perpendicular to the plane of symmetry). 
If the character of the one- dimensional irreducible representation is 1 (-1) the boundary 
condition is 7/(0) = (7/(0) = 0). When the integral equations for two components of the 
order parameter decouple, the boundary conditions can be similarly written. 

The irreducible representations and corresponding basis functions of are listed in 
Table I. We assume that the spins in the triplet state are ordered along z axis, this is 
denoted by the spin matrix c. In this particular case the equations for the order parameter 
can be reduced to Eqs.(10) as well. 

When the angle 9q between the normal to the boundary n and z axis constitutes tt/2 
and the azimuthal angle 0o (between x and n) equals fc7r/6 (k is integer) the boundary 
conditions for different representations are enlisted in table II. The notations 

f 7/i = 7/1 cos 0o + 7/2 sin 0o , 12 \ 

\ 7/ 2 = -7/i Sin O + 7/2 COS O 

are introduced there for E\ g , E\ u representations. 
Similarly for E 2g , E 2u representations: 

f t)i = 7/i cos(20 o ) - 7/2 sin(20 o ) 
\ 7)2 = T)i sin(20 o ) + 7/2 cos(20 o ) 

If n coincides with z we obtain the following boundary conditions for one-dimensional 
representations and for both components of two-dimensional representations : for A lu , A 2u , 
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Big, B 2g , Ei g , E 2u - rj(0) = 0; for the remaining representations - rf(0) = 0. We see that 
the microscopic structure of the gap function reveals itself in boundary conditions to the 
Ginzburg-Landau equations, although these equations may be similar for different represen- 
tations. 

Now let's discuss the solution of Eq.(10) for arbitrary orientations. It is easy to check, 
that for the basis functions of Table I in the case of spherical Fermi-surface these integral 
equations decouple, if transformations flI2D,flI3p are applied for E\, E 2 functions respectively: 

m = fi ^ 2C E f / d 2 smup){^eM-\—\\i-i'\)+ 

J 1pt{p)d 2 S m J ° J N V l 

+ *tMeM-\—\(l + l'))}dl'. (14) 

M Vl 

There is no summation over subscript i in this equations and they look exactly like equations 
for the one-component order parameter. The functions ipi(p) are related to ipi(p) by the same 



equations (|12|),(|13|). The gap function may be rewritten as A(p, I) = ipi(p)fji(l) +^2(p)v^) 
At the distances £ ^ I <C £{T) from the boundary Eq.(ffl]) has a solution fji(l) = 
C(l + qi), the value of may be obtained from the variational procedure developed in 
which usually turns out to be a good approximation. The exact relation for is: 

a>(:1, ' / ' A f tf(P)v?d 2 S = £=( Mp)mP)+HP))h\ 3 d 2 S+[D mm [ tf{p)d 2 S]-\ 



7T 2 T 2 J Tt ^' 1 48T 3 

(15) 

where D m i n is the minimal value of the functional: 

(™ q{l)[q{l) - ^ K{U)q{l')dV]dl 

[Jo°° i(i) E ( 1 )} 2 ' 1 ' 

Here K(l,l') is the kernel of the original integral equation flTjj ) (fji(l) = / °° K(l,l')rji(l')dl') 
and 

EVA T\V^ f 7 f-\$i(P) +^i(P))KI / \ 2uJ ™\l\j2o f,n\ 

E(l) =vrTA> / ipi(p) - exp(-| \l)d S. (17) 

The function q{l) in Eq . (|16D is supposed to tend to the constant value, when I — > oo. 
Estimating (|16|) for constant q(l) we obtain approximate values of g«: 



( 71 
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\336((3)T J Vl >o 



\Mp) + Mp)Yv?d 2 s+ 

+ (18) 

4n T J Vi >Mp) -A{p)] 2 vid 2 S J J 

where the subscript v \ > denotes that the integration is made over that part of the Fermi 
surface, where Vi > 0. This equation may be applied for one-dimensional representations 
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directly for the boundary condition qr}'(0) = T](0) with the non-til ded values 77, since no de- 
coupling problem arises in this case. Eq. fll8|) remains to be valid for complex basis functions 
as well after simple replacement of squares of the basis functions and of sums and differences 
of the basis functions by squares of the corresponding moduli. 

Particular values of qi, generally speaking, essentially depend upon the explicit particular 
form of basis functions ipi(p), which may differ even if they belong to the same symmetry 
representation. Only values q = 0, 00 governed by the symmetry are invariant for the given 
type of pairing (irrespective of the shape of the Fermi-surface of the given symmetry). For 
the specific basis functions of E lg type from Table I the integration in flTg| ) results in: 



,s ( 95tt 3 , 12C(3) 



q{ = — sin 2 29 — + % ' tan z 29 

Hl T c °V18432C(3) 5vr 3 

_ v F sin 2 6 / 5vr 3 4C(3) \ 
52 - % 1 + 2cos^ l^36CM + ~5^ tm H (19) 

The boundary condition takes the form ?i^(0) = fji(0). So if 9 = 7r/2,</> = we obtain 
77i(0) = (since q\ vanishes), while rj' 2 (0) = (since q 2 = 00). Similarly other orientations 
from Table II and the case of n along z may be checked. From (19) it follows that beyond the 
relatively narrow vicinities ~ (^o/^iT)) 1 ^ 2 of the directions for which qi = 00, the parameter 
qi has the order of £ , because near these directions t& ~ £(0)/(A# ) 2 , where A9 is the 
angular displacement from these directions. 

For the triplet pairing of .Ei-type Eq.flPSD and Table I yield: 



t Vf cos 2 9 ( 5tt 3 14C(3) 2 
qi T c 3 sin 2 #0 + cos 2 #0 V 1344C(3) + 15vr 3 C 

^ = 00. (20) 

If the functions are those from Table I of E 2 type we obtain in the case of singlet pairing: 

5 2 9 f 5tt 3 
T c 1 + 2 sin 2 9 U536C(3) 



_ v F cos 2 9 / 5tt 3 4C(3) 2 \ 
qi ~~^ + 2 sin 2 #A 153603) + ^r r H 



vf ( 5tt 3 , „„ . . 7C(3)(3sin 4 ^ -sin 2 ^ + l) 2 

T c (l + 2 sin 9 ) \6144(,(3) hn 6 sin 29 

(21) 

And in the case of triplet pairing: 

t v F 2 „ ( 19tt 3 32C(3) „ „ 

* = t M + 6 ^C- wa. (loAsy' 100 - 276si " 29 » + 195si " 49 » )+ 

32C(3) (22 - 59 sin 2 9 + 41 sin 4 9 Q f tan 2 9 \ 
+ 15tt 3 4- 12 sin 2 fl + 27sin 4 9 )' ^ ' 
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The values of qi enable to write down the surface free energy and then (considering 
the total free energy) to obtain the general boundary conditions for the Ginzburg-Landau 
equations, not necessarily for the case, when the order parameter is changing only in the 
direction normal to the boundary. In the weak coupling limit the gradient terms of the 
Ginzburg-Landau Hamiltonian, as is known, take the form 



So in the case of spherical Fermi-surface and functions from Table I we get, that for the 
^i-type of pairing K\ = K 2 = K 3 = K, besides in the case of singlet pairing K 4 = 3K, 
while for triplet K± = K. Taking into account these relations one can find the explicit 
expressions for hi n ,h2 n in (3). Then comparing Eq.(3) with the microscopic boundary 
conditions ^^'(0) = 77^(0) one obtains the coefficients 6,j resulting in the following expressions 
for the surface free energy: 



3K . ± . ± , a K(sin 2 flo + 3cos 2 #o), ± . ± 

F surf = ~rTa~\\ r ) 1 cos «Po + V2 sm O I + jt^t \V2 cos 0o - 771 sm 1 , 
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, K (3 sin 2 8 + cos 2 6» ) . . , „ 

F surf = — tJTn ^Ir/icos^o + ^sm^l 2 . (24) 



Here g*'* are given by (|I~9l) , (|20l) respectively. Note, that free energy is written down in terms 
of original 771, 772. 

If the pairing is of E 2 type, then for the spherical Fermi-surface one has K 2 = K 3 = 
and K± = K/3 (K4 = K) for singlet (triplet) pairing. Thus the surface free energy acquires 
the form: 

F s surf = K ("sin 2 6»o + C ° S Q — ) ( s ] . | i]i cos 2(j) -r] 2 sin 20 o 1 2 + S ] Q \ I V2 cos 20 o +Vi sin 20 o 1 2 } , 
\ o / l qi{t>o ) <l2\yo) > 

F\ urf = k{ 1 |r/icos20o - r/ 2 sin 20 o | 2 + t ] |r/ 2 cos 20 o + r]i sin 20 o | 2 |. (25) 

Like the specific values of the parameter t&, the particular form of the orientational 
dependence of the surface free energy is governed by the explicit particular form of the basis 
functions. Any angular dependence of qi and F sur f can be certainly expressed in terms of 
n 2 , i.e. sin 2 6q = n 2 + n 2 . The surface free energy can be represented as a simple quadratic 
form over n 2 ,n 2 ,nl for a special kind of basis functions only. According to (24), (25) higher 
powers of n 2 are also involved in the expression for F sur f for the basis functions from Table 
I. If the free energy (1) is invariant under the rotations on arbitrary angles around z-axis 
and qi depend only upon 9 , one can write down volume and surface free energies in a O - 
independent form simultaneously. It takes place for a particular choice of basis functions 
for two-dimensional representations indicated in Table I and for rotationally invariant Fermi 
surface. So, in the particular cases (1), (24) or (25) the boundary conditions in fact do 
not result in a o -d e P en dence, since fji can be used in the total free energy. One can easily 
demonstrate that for both two-dimensional representations there are basis functions breaking 
rotational symmetry of the total free energy functional. 
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If the value of qi is of the order of £o, we have fji(0) ~ fjioo (£,o/£,(T)) <C fj ioo at the 
boundary, where fj ioo is the value in the depth of the superconductor. So near T c we can set 
approximately fji(0) = 0. Hence, as it follows from (19), (21), (22), the prevailing boundary 
condition (i.e. for the most of orientations of the boundary) for E^-type of symmetry (both 
triplet and singlet) and for singlet pairing of E\- type is: 

r?i(0) = 77 2 (0) = 0. (26) 

For the specific triplet basis functions of Ei-type from Table I the prevailing boundary 
condition is: 

?7i(0) = ^(0) = 0, (27) 

but if other basis functions without rotational symmetry are used, or non-spherical Fermi- 
surface is considered, the prevailing boundary condition will be Eq. (p6|) for this case as 
well. 

If the scattering at the boundary is diffusive, rather than specular, the appropriate 
values of will be of the order of £ as is shown in P9"| , so we can use Eq.(^) as a 
good approximation near T c for all orientations in this case (and for all anisotropic pairings 
irrespective of the crystalline symmetry). 



III. TWO-COMPONENT ORDER PARAMETER: H C 2{T) FOR A 
SUPERCONDUCTING FILM 

Let the crystal axis of high symmetry (axis z) be parallel to the boundary of a film 
with thickness 2d, axis x be perpendicular to the film surface and a magnetic field H \\ x 
(see fig.l). Though our goal is to describe some properties of a hexagonal superconductor 
UPts as one of a probable candidate for a two-component unconventional superconductor, 
there are some formal reasons (associated with the simplifying of the problem) to consider 
firstly a more general case of a tetragonal superconductor. Besides the case of a tetragonal 
unconventional superconductor is of interest as itself, since the question about the type 
of pairing remains unsolved also in some tetragonal heavy fermion superconductors like 
CeCu2Si2 and URU2S12. Then the Ginzburg-Landau free energy up to the second order 
terms reads 

d 

F = J J dydzl J dx [a^r]* + Kxp^p^ + K 2 p*r]tPjVj + K 3 p*VjPjVi + K *P*zV*PzVi+ 

~d 

+K 5 (\p x r] 2 \ 2 + |2V?i| 2 )] +Si (M 2 |,=_ d +M 2 \ x=d ) +S 2 (M 2 \ x =-d +M 2 U=d)} (28) 

The terms describing the surface contribution to the free energy in Eq.(^), correspond to 
taking into account the boundary conditions for the order parameter. We have not written 
one more surface term S^r]^ + ^^2) since for the given orientation of crystal axes relative 
to the film boundary one has £3 = 0. The coefficients Si, £2, S3 may be determined by 
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comparing the boundary conditions following from (|28|) with those based on a microscopic 
theory. 

It is convenient to choose the vector potential in the form A z = Hy and to introduce 
dimensionless quantities e = (K 2 + K 3 )/K, \i = (K 2 — K 3 )/K, ( = 2K$/K, v = (K 2 + 
K 3 - K 5 )/K, a = -ac( 1/2 /(2hK 4 \e\H), S 1>2 = 2S^ 2 l/{Ti 2 K), x -> xl, y -> yl, z -> zl, where 

(hc/2\e\Hy/ 2 . Then Eq . fl2"8|) is transformed as follows 



K = 2K X + K 2 + K 3 + K 5 , I = C 

L 



-1/4 



-KhH 

2 



dydzl / c/x[-5(|r/ 1 | 2 + |r ?2 | 2 ) + (l + z/)(|^ 1 | 2 + |^ 2 | 2 ) + 



+(1 - v){\d x 7] 2 \ 2 + \9yVi\ 2 ) + y 2 (\Vi\ 2 + \V2\ 2 ) + (e + n){d x r,ld y r] 2 + c.c.)+ 

+ (e - v)(d x r]*d y r}i + c.c.)] + St (jr]i\ 2 \ x= - L + \vi\ 2 \x=l) + S 2 (\r] 2 \ 2 \ x= ^ L + M 2 |*=l)} (29) 

Here 2L = 2d /I is a dimensionless film thickness. One should take into account that for the 
given gauge of A it is possible to choose the order parameter 77 = (7/1,7/2) to be independent 
of z in considering the upper critical field (it corresponds to the choice of the origin at the 
y-axis) . 

Variation of the free energy functional ( p9f ) on 77*, rj 2 results in the equations 
-{l + v)dl + y 2 -(l-v)df 



7/1 - 2ed xy r] 2 = arji 



[1 - u)d 2 + y 2 - (1 + v)d 2 7/2 - 2ed xy r} x = ar] 2 



(30) 



and in the boundary conditions 

1 + v)d x r)i + (e + [i)d y r) 2 ± ^t/i 

1 - v)d x r} 2 + {e- f^dyij! ± S2V2 



\x=±L 



\x=±L 



(31) 







Due to linearity of Eqs.([3(], ^T]) one can consider only those solutions which are charac- 
terized by the definite values of parities relative to the independent inversions of x and y 
axes. Then it is easy to see that the components rji and 7/2 belonging to the same solution 
of Eqs.([3(], [Tl]), have opposite values of parities both on x and on y. 

The condition of positive definiteness of the sum of the gradient terms in (|29|) results in 
the following restrictions 



1 ± v > \e ± fx\ 
In particular, it follows from (R3) that |e| < 1 



C>0 



(32) 



\fi\ < 1 , \u\ < 1. Thus one can usually 



regard the parameters e, /x, v as quite small (for the majority of their admissible values). We 
shall employ below the perturbation theory up to the second order in e for the solution of the 
problem of the upper critical field for the film. For a hexagonal crystal one has K 5 = and 
hence v — e. However, for the simplest qualitative examining of the problem it is convenient 
to assume the parameters v and e to be independent and consider only e as a small quantity. 
It is, generally speaking, possible for a tetragonal crystal. The problem of determining of 
the minimal eigenvalue a will be solved for specular interfaces and for diffusive boundaries 
separately. 
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A. Specularly reflecting film surfaces 



Since yz-plane is a symmetry plane both for the symmetry group and for D^, for 
a given crystal orientations and two-dimensional representations of those groups specularly 
reflecting film boundaries correspond to the limit Si — ► +00, S% = in a phenomenological 
Eq. (pi)) . Then the boundary conditions read 

Vl\x=±L = , d x r] 2 \ x= ±L = (33) 
The solutions of Eq.(^) with boundary conditions ( p3|) can be represented as follows 

' 7] ln = COs(k n X + ipn)r]ln{y) 

n = , ±1 , ±2 , • • • (34) 

^ r] 2n = Sm(k n X + (p n )V2n(y) 

where k n = (n/2L)n, ip n = (7r/2)(n + l) and functions r]i n (y) , r)2 n {y) must satisfy the system 
of equations 

(y 2 - (1 - u)d 2 )rj ln (y) - 2ek n d y rj 2n (y) = (a - k 2 n (l + v))-q ln (v) 

n^O (35) 

(y 2 - (1 + v)d 2 y )r] 2n {y) + 2ek n d y rj ln (y) = (h - k 2 n (l - iy))rj 2n (y) 

(y 2 -(l + iy)d 2 ) V 2 (y) = ~ar ]20 (y) , n = (36) 

The problem of the upper critical field becomes especially simple and transparent in the 
zero order approximation in the parameter e, when Eqs.(^) for r]x n , r]2 n decouple. From the 
right sides of Eqs. Q55]) one can see, that dependence of any component of the order parameter 
upon x results in decreasing of corresponding "critical temperature" T n (as compared to T c in 
the absence of a field). As a result all the quantities H n (T), following from fl35| ) for n 7^ as 
candidates for the upper critical field, describe the magnetic fields decreasing with increasing 
temperature and vanishing at T n < T c . Furthermore, it follows from the equations for r]2 n 
that the maximum value of magnetic field is reached for the equation ( p6|) with n — 0, since 
a, k 2 oc 1/H and v < 1. As T = T c > T n {n 7^ 0) the corresponding magnetic field H Q {T) 
must coincide with H c2 close enough to T c . Then one has for the slope of the upper critical 
field H' c2 oc 1/(1 + v) 1 / 2 . At the same time it follows from the equations for r] ln (y) that 
the slopes of corresponding temperature dependent magnetic fields oc 1/(1 — u) l l 2 . Hence 
for v > these magnetic fields increase with decreasing of temperature faster than the field 
Hq(T), defined by (|36|). Since the equation with n = 1 leads to the maximum admissible 
value of the magnetic field among all the equations for r] ln (n ^ 0), one comes eventually to 
the following result at e = 0. In the case v < and for the specularly reflecting boundaries 
the upper critical field for the film coincides with the one for the massive sample for all 
temperatures and is described according to (j36|) : 



H c2 {T) = % (37) 

M 2 (T)J<(1 + v) 
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where £ 2 (T) = h 2 K/2\a\ = - T/T c ) , $ = nhc/\e\ . 

But as far the case v > is concerned the expression (B7) describes the upper critical 



field only within the certain temperature interval T k 
H C 2(T) suffers a kink and is described as follows 



(0) 



<T <T C , and for T <T k 



(0) 



^ 2 0) (T) 



2<oVC(l 



1 



T 

2; 



Ad 2 



;i + 1/) 



the curve 



(38) 



For the point of kink one has from (|37|) 



(0) 



(1 + z/) 3 / 2 



4d 2 (x/TT^ - v 7 !" 



z/ 



(39) 



In order to satisfy the condition {T^ — T c \ <C T c the inequality d 2 v 3> £ 2 must hold. 

Since the solutions with n = and n = 1 have a different symmetry- opposite parities 
relative to the inversions on x and on y- for any value of e (and not only for e = 0), it 
is obvious that the kink at the curve H C 2(T) must remain in the case e ^ as well ( for 
the given orientations of magnetic field and crystal axes). In the case e ^ the quanti- 
tative consideration of the problem becomes more complicated. However, as for hexagonal 
superconductors v — e, it should be carried out in order to describe the effects in question 
correctly. Corresponding calculations are represented in Appendix. According to the solu- 
tion of Eq.(|3~5"D for e ^ 0, the slope H[(T) becomes temperature dependent even within the 
region of applicability of the Ginzburg-Landau theory T c — T T c . At T k < T < T c the 



slope coincides with the bulk result fl37|) . At T < T k one gets in accordance with (A. 25) up 
to the second order terms in the parameter e inclusively 



T, 



$ 



1 + v- 



X 1/2 d 2 

e 2 {l-C 2 f' 2 } 



dt 



(l-„2)l/2 J (l-CH 2 ) 3 / 2 



t f3(»)-j(u,H) 



(40) 



where 



C 



y/l + U- v 7 ! - 



VI + V + VI - V 



'1-1/ 
1 + 1/ 



7(1/, 10 



7TZ/$r 



8C 1 / 2 rf 2 (l + z/) 1 /2// 



(41) 



Under the condition e 2 , |i/| < L 2 = 2ird 2 ( 1 / 2 Hi(T) /<J> the upper critical field below the 
kink temperature takes the form 



# c2 (T) 



2< 2 C 1 / 2 ^ 



7/ 



T„ 4rf 2 



2 ( 2 £ 2^_ £2)3/2 j dyy m 



(1-1,2)1/2 J ( 1 _ C 2 y 2)3/2 







(42) 



In the limit e 2 /L 2 — > 0, z//L 2 — > the expression (pE2]), becomes exact. 
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The crossing point of the right lines (0), (0) is 



^-^j 1 -4d T ^T^_^T^[ i + ,y (i-^)va J (i_ (7V) 3/2 , 



(43) 



It is assumed in ( pf) that the kink point (Tk,Hk) lies in the region where the curve 
H\(T) nearly coincides with its asymptotic straight line. It is valid if ^qE 2 / (2Txd 2 C > 1 / 2 ) <C 
which in fact is equivalent to the condition ve 2 <C 1. For e — v the Eqs. (0), (|42|), (|43|) 
describe the upper critical field in the film for a hexagonal unconventional superconductor 
with two-component order parameter. The corresponding dependence of H c2 (T) is presented 
in fig. 2 (line 1). 

Considering the upper critical field for some other orientations of crystalline axes and 
then comparing the result with one described just now, one can get several simple tests for 
identification of one or another type of pairing. For instance, one can distinguish between 
the order parameters transforming according two-dimensional representations E\ g , E\ u and 
E 2g ,E 2u , by examining the temperature dependence of the upper critical field near T c in 
the film with specularly reflecting boundaries for two different orientations. Since for a 
hexagonal crystal the x- and y-directions within the basal plane, generally speaking, are not 
equivalent, let now axis y be perpendicular to the film surfaces, axis x be parallel to the 
boundary and a magnetic field H || y. The results for this particular case can be readily 
obtained from the expressions fl3~7|), ([40]), if one takes into account that initial Eqs.(|30|), 



(j33|) must be transformed to the case along the following way. In the Ginzburg-Landau 
equations ([30]), written in the crystalline coordinate system, the only modification to the 
present case is the substitution y 2 — > x 2 , which corresponds to the field direction along the 



y crystalline axis. However, as far as the boundary conditions (33) are concerned, different 
modifications should be used for the representations E\ g , E\ u and E 2g , E 2u . In the case 
of the representations E\ g ,E\ u the order parameter components, according to Eqs.(|i~2|) at 
<po = 7r/2, must be interchanged in the boundary conditions (^) along with d x — ► d y . Then 
the resulting equations combined with the boundary conditions exactly coincide with the 
initial ones after the formal substitution x <-> y, r\\ <-> 1] 2 . Hence the upper critical field 
remains the same for both orientations in the case of E± g , E\ u - representations, but it is 
not the case for the representations E 2g ,E 2u - Indeed, according to Eqs.flHID at 0o = tt/2 in 
the latter case the boundary conditions take the same form (^) but for the substitution 
dx d y . Then after the formal substitution x — > y, y — > —x the resulting equations differ 
from the initial equations only by the sign of e and v (for a hexagonal crystal v = e ). Since 
as was shown above, for specularly reflecting boundaries the kink occurs only in the case 
v > 0, we are coming to the assertion that if the order parameter transforms according to 
the representations E 2g , E 2u , the kink must be present for one of the orientations and be 
absent for the other. This statement holds also for <f) = (2k + 1)tt/6, where k = 0,1, 2, .... 

For the specific case of a tetragonal superconductor with two-component order parameter 
one can similarly suggest the following test permitting to distinguish relative signs of the 
parameters e and v. Indeed, let the normal to the boundary n lie between the crystalline 
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axes x,y constituting angle 7r/4 with each of them, and H || n. The film boundary coincides 
in this case with one of the symmetry planes of the group D^. One can show, that the 
Ginzburg-Landau equations and the boundary conditions for this particular crystal orienta- 
tion are exactly the same as Eqs.(|30D, (|33D after the substitution v e, if they are written 
for the functions r\\ ± r? 2 in the variables (x ± y)j y/2. Then, if the parameters v and e have 
opposite signs, the kink in the temperature dependence of H c2 (T) occurs only for one of the 
orientations. In the other case, when the parameters are of the same sign the kink is present 
(absent) for both orientations simultaneously for v > {y < 0). 

Let now temperatures T cl and T c2 in coefficients Oj = a(T — T c i ), i — 1, 2 be not equal 
to each other and AT C = |T cl — T c2 | <C T&. Then one has instead of a(|r/i| 2 + \i]2\ 2 ) the 
terms cii|^i| 2 + Q-2 1?72 1 2 in the expression (p8|) for the free energy. Similar to consideration 
given above one can get that at T cl < T c2 and v > the kink of the curve H&iT) occurs 
both in a massive homogeneous sample and in a film. The boundary conditions significantly 
influence the kink location in the case £q/gP ^ AT C /T C . On the other hand at T c \ > T c2 and 
v > the kink doesn't exist in a massive homogeneous sample but appears in a sufficiently 
thin film when Ti(d) < T c2 . Furthermore in the case v < and T cl < T c2 there is no kink of 
H C 2(T) both in a massive homogeneous sample and in a thin film. At last in the case v < 
and T c i > T c2 there is a kink in a massive homogeneous sample while it is not present in the 
case of sufficiently thin film (at Ti(d) < T c2 ). 



B. Diffusive electron scattering from the boundaries 

While under the condition of a specular reflection of electrons the kink of if c2 (T) exists 
only in the case v > 0, it will be shown below that in the case of diffusively scattering film 
boundaries the kink of if c2 (T) always takes place no matter what is a value (and a sign) of 
the parameter v. This assertion can be easily understood qualitatively from the following 
simple consideration. 

According to the results of Section II, in the case of diffusively scattering film boundaries 
the boundary conditions read 

Vi\x=±l = , r) 2 \ x=±L = (44) 

In zero-order approximation in a small parameter e the equations (|30|) decouple and 
their solutions with account of the boundary conditions ([0]) can be represented in the form 
T)in(x, y) = cos(fc„x + ri )?7 l - ( ^(?/). Solutions for rj^(x,y) may be obtained from the solutions 

T)i^(x,y) by changing v — > —v. In contrast to ( |33D the boundary conditions ( 44]) don't allow 
the value n = to be admissible both for r)i n and for r] 2n . Then the least changing of a 
critical temperature due to the boundary conditions is realized for n — 1 and is proportional 
to a corresponding coefficient in front of the term with the second order derivative over the 
coordinate x in Eq.(^), so that 

T c -T x l + v 



T r -To \-v 



(45) 



On the other hand one can see from the equations for tji and r]^ that a ratio of slope 
H{ 0) '(T) (for 7]{ 0) ) to a slope Hf Y {T) (for r/f) is 
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H[ oy (T) _ (l + v\V 2 



v 



(46) 



Thus at v > the field H 2 (T) drops to zero at temperature T 2 greater than T\ which 
corresponds in turn to zero value of the field H\{T\). However, since a slope H\ (T) has 
a greater steepness there appears a crossing of curves H[ 0) (T) and H { 2 0) {T) which results 
in the kink of H c2 (T). Under the condition v < one gets from (§5|), fl4"6|) T 2 < Ti while 



1-^2°^ CO I > \Hi V> CO I- It leads again to the crossing of the curves H^^T) and H^'iT), 
and hence to the appearing of the kink. Thus in zero-order approximation in the parameter 
e it follows from Eqs. @, (@) 



r(o)' 



(0), 



r(0). 



^i 2 0) (T) 



2tt^C 1/2 [l + k|^n(T-T ? 



(oh 



1/2 



T 

l -T, 



Ad 2 



[l-\v\sgn(T-T^)) 



(47) 



(0) 



4d : 



f(2 + v / T^) 



(48) 



H, 



(0) 



7T$ r 



8d 2 C 1/2 



V 



As it follows from Appendix (see (A. 31)), for e ^ the temperature dependence of the 
upper critical field is described up to the second order terms in e inclusively by the equation 



T 2ttC 1/2 £ 2 / x1/ o tt 2 £ 2 / 



V- 



4e 2 (l-C 2 ) 3 / 2 } 



7T 



'1 - Z/ 2 )V2 



dt- 



(49) 



4fc 2 ^f - l) /(16C 1/2 rf 2 //). 



where z/ = |i/|s<?n(T — Tfc), <5(fc, z/) = 7r$o\/l — v 

Thus at e ^ and diffusively scattering boundaries the upper critical field H c2 (T) ex- 
hibits a feebly marked nonlinear temperature dependence (in measure of the parameter e 2 
usually multiplied by a small numerical factor) for any sign of v. Under the condition 
1 < L 2 = 2ix( > l l 2 d 2 H/Qo one can get the exact asymptotic behaviour for the upper critical 
field 



H c2 {T) 




2vrCn 2 (l 



Ad 2 



e 2 (l - C 2 f' 2 

(1 - v 2 yi 2 



dt- 



t P(-o) 



CH 2 fl 2 



(50) 
obtained 



For T <T k Eq . (|50|) completely coincides with the asymptotic expression ([ 
for the case of specularly reflecting boundaries. 

In the opposite limit 1 ^> L 2 one has up to the second order terms in e inclusively (see 
(A.37)) 
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„ m go tJjzjtSlzfh IVI) 

- l ' - 2, C "^(i + w (i - ^ (i - ^Foot ( f ys)- 



In Eq.(pTf), which is applicable just below T c , may be changed on v since it follows 



from (|48|) that the kink point lies at L ~ 1. 

The temperature dependence of H c2 (T), described by (49) is presented in fig. 2 (line 2). 

If T cl ^ T c2 then, like in the case of specularly reflecting boundaries, one can easily com- 
pare the conditions for the presence or absence of a kink of H c2 (T) in a massive homogeneous 
sample and in a thin film with diffusively scattering boundaries. A distinction from the case 
of specular reflection appears here at v < 0. Then for T c \ < T c2 the kink is absent in a 
massive homogeneous sample, but takes place in a sufficiently thin film with Ti(d) > T 2 (d). 
For T cl > T c2 and v < a kink occurs for the both cases. 



IV. H C2 (T) FOR A SUPERCONDUCTING FILM IN THE CASE OF 
ACCIDENTAL DEGENERACY 

In the case of accidental degeneracy two one-component superconducting order parame- 
ters with different symmetry properties correspond to critical temperatures T cl , T c2 , which 
are very close to each other (|T cl — T c2 | = AT C <C T cl ). The Ginzburg-Landau functional 
for uniaxial superconductor in such a situation may be represented up to the second order 
invariants as follows 

r d 

F = J dydz I J dx ciiM 2 + a 2 \r) 2 \ 2 + K^p^ 2 + |jy?i| 2 ) + K '(\p x rj 2 \ 2 + \p y r)2\ 2 ) + 
{ d 

+K 4 \p z r ]l \ 2 + K' 4 \p z r] 2 \ 2 ] + S 1 (Mf x= _ d + MjUd) + S 2 (l^ljU-d + MjL=d)} (52) 

Here a\ = a± (T — T cl ) , a 2 = a 2 (T — T c2 ). The coherence lengths ^ = h 2 K/aiT c \, 
^ ± = h 2 K 4 / ai T cl} #2 = h 2 K'/a 2 T c2) ^o 2 ± = h 2 K'Ja 2 T c2 will be used below. 

In contrast to (^), where there are gradient terms containing d x r]id y 7]j with i ^ j, the 
terms of such kind are absent in (|52|) . In consequence of this fact two Ginzburg-Landau 
equations for rji and 7] 2 decouple, if the terms of the fourth order are neglected in the 
expression for the free energy . This circumstance essentially simplifies the problem of the 
upper critical field in a film. Let's consider the same as in a previous section orientations 
of the crystalline axes and of the magnetic field relative to the film boundaries, and the 
same gauge for the vector-potential A z = Hy. Let the boundary conditions for the order 
parameters r\\ , r\ 2 be [fH]] 

Vi\x=±d = , d x T) 2 \x=±d = (53) 

According to the results of Sec. II, these boundary conditions are realized for the given ori- 
entations, in particular, for specularly reflecting boundaries of a hexagonal superconductor, 
if the order parameter rji corresponds to one of the one- dimensional irreducible representa- 
tions A 2g , A 2u , B 2g , B 2u of D^h and the order parameter r\ 2 - to one of the representations 
A-xg-, Ai u , Big, B\ u . 
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The upper critical field may be represented as H c2 = max(ifi(T), H 2 (T)), where the 
linear functions of temperature Hi(T) and H 2 (T) are associated with the solutions of the 
Ginzburg-Landau equations rji = 0, r\ 2 = 7/20(7/) and 971 = cos(7nc/2<f)7/i (7/), r/ 2 = re- 
spectively. Here 7/10(7/) and 7720(2/) are related to the zero Landau levels of corresponding 
equations, and Hi(T),H 2 (T) - the appropriate magnetic fields. 

Firstly neglecting the difference between T c \ and T c2 , one has T 2 = T c and ATi = 
T c - Ti oc Then H[{T)/H' 2 {T) = £ ±£o|| /£<u£o|| • Thus the inclined straight line #i(T) 
comes to zero at a lower temperature than H 2 (T) and with decreasing of the temperature 
crosses the line H 2 (T) only under the condition £ ||£ _|_ > £o||£o_i_- If T cl < T c2 , the above 
condition has the form ^C' ±T c2 > £o||£o±Tci and remains for the massive sample as well. 
The influence of the film boundaries on the crossing point position becomes noticeable if 
£o||/^ 2 ~ AT C /T C . For T c i > T c2 the kink of H c2 (T) is absent for a massive sample and appears 
for a thin film if Ti(d) < T c2 . Under the condition £ ||£o±^c2 < £o||£o_L^ci and T cl > T c2 the 
kink occurs for a massive sample and dissapears with decreasing a film thickness, when the 
condition T\{d) < T c2 becomes valid. 

In the case of diffusively scattering boundaries the boundary conditions are 

T)i\ x =±d = , r]2\x=±d = 0. (54) 

As just above, let's firstly neglect the difference between T c i, T C 2 and compare linear 
functions ifi(T), H 2 {T) which are related with the solutions 7/1 = cos (nx/2d) 7710(7/) , r\ 2 = 
and r/i = , r/ 2 = cos (7tx/2c/,) 7/20(7/) respectively. These solutions satisfy the Ginzburg- 
Landau equations following from ( |52D together with the boundary conditions (|54]). Here 
H c2 = max(ifi(T), H 2 {T)). Since for the case in question one has 

^Ei = mi = Mi / 55 n 

the condition for H c2 (T) to have a kink becomes much more restrictive than in the case of 
specular reflection. Indeed, if £ || > £o||> ^ ne kink takes place under the condition £ 0± /£o_i_ > 
fo||/fo|| > L In the opposite case ^ || < Co|| the kink occurs if £ ±/£o_l < Co||/Cd|| < L The 
latter condition means, in particular, that if the relation £ 0± < holds for the coherence 
lengths in the perpendicular to the basal plane direction, the opposite relation £L > £ || 
must be valid for the coherence lengths in the parallel direction. Though the relationship of 
such kind is, in principle, admissible, it is a quite strict demand MM. It is in contrast with 



the result of previous section for two-component order parameter and the same boundary 
conditions (|54D , that the kink in the temperature dependence of the upper critical field 
always takes place for a thin film. 

If T ci T c2 and, for example, T c i > T c2 , then the kink occuring for a massive sample 
under the condition CLC'o±T c2 < £o||£o_lT c i, remains at £o||-^ cl < £oj|-^ c2 f° r a thin film with 
diffusively scattering boundaries too. In the opposite case £o|l^ci > £o||^c2 the kink will 
disappear with decreasing of film thickness, when T2(<i) > T\{d). 
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V. CONCLUSIONS 



Boundary conditions for the unconventional superconducting order parameter have been 
derived on the basis of microscopic theory for specularly reflecting boundaries. Using these 
boundary conditions (and those for diffusive scattering) for the problem of the upper critical 
field in a superconducting film, the appearance of a kink in a temperature dependence of 
H c2 (T) has been examined. The presence or absence of a kink depends upon the ratios of 
coefficients at gradient terms in the Ginzburg-Landau functional and upon the quality of 
boundaries (specular reflection or diffusive scattering). If due to some other reasons the 
kink in a temperature dependence of H c2 (T) takes place even for a massive homogeneous 
superconducting sample, the study of the location of a kink point in a thin film from the 
same compound may give a useful information about the type of pairing and the relations 
between the Ginzburg-Landau coefficients. 

In the case of diffusively scattering boundaries and a two-component superconducting 
order parameter the kink turns out to be always present for the crystalline axes and magnetic 
field orientations considered above. Then the location of a kink point is qualitatively defined 
by the condition d ~ £(T). At the same time in the case of accidental degeneracy the kink 
of H&iT) in a thin film occurs due to the influence of diffusively scattering boundaries only 
under quite hard restrictions for the relative values of the coefficients at gradient terms. 

For a two-component order parameter and e 7^ the kink, strictly speaking, takes place 
only for considered above orientation of a magnetic field perpendicular to a crystalline axis 
of high symmetry. But for small enough values of e and \i a smooth behaviour of the curve 
H C 2(T) may exhibit an abrupt change of slope resembling the kink very much. 

Recently incommensurate structural macroscopic modulations in heavy fermion UPt^ 
have been observed along several crystalline directions |29 ]. Proceeding from the consider- 
ation of a thin film as the simplest model for studying the influence of an inhomogeneity 
on the sample properties, it is natural to suppose that for the appropriate orientations of a 
magnetic field smooth inhomogeneities of macroscopic scale may lead to the kink of H&(T) 
as well. Note, that in the case of a thin superconducting film and of two-component order 
parameter the splitting of the superconducting transition in the absence of magnetic field 
may also be provoked due to the boundary conditions [l],[42|,[|3|] . So, it is possible to consider 
the structural macroscopic modulations as one of the probable sources for the kink of H C 2(T) 
and for the splitting of the critical temperature in UPtz- 

After this article had already been submitted a paper by K.V.Samokhin appeared in 
Zh.Eksp.Teor.Fiz. 107, 906 (1995), where the boundary conditions for unconventional su- 
perconductors were also considered and applied to the problem of the surface critical field 

Hc3- 

The research described in this publication was made possible in part by Grant No 
M2D000 from the International Science Foundation. 



APPENDIX 

It is convenient to look for the minimal eigenvalue a of the problem (^) using the 
perturbation theory in the parameter e and transforming the problem to the form 
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H V + = a* (Al) 

where 

v=( ~ 2£d *A 

' {m)' \-2ed xy )> 

-(l-^ + y 2 -(l + ^ 2 J- (A2) 

In terms of the operators 

^ _ iy/1 - vd y + ^ + _ VI - ^ - iy ^ _ iy/1 + vd y + iy ^ + _ iy/1 + vd y - iy 

y/2(l-v) 1 / 4 ' " V2(l-v) 1/4 ' ~ v / 2(l + ^) 1/4 ' " v / 2(l + ^) 1/4 

(A3) 

which satisfy the conditions [A, A + ] = 1, [B,B + ] = 1, the unperturbed Hamiltonian reads 

fr _ ( "(I + ^ + 2vO^(i+i + \) \ 

° ^ -(1 - v)dl + 2y/T+l}(B + B + \)) 

From Eq.(A.4) one can see that a complete orthonormal set of eigenf unctions of the 
operator H may be represented as a set {^ mn } = {^mm ^mn}-> where 

*l=(f") , ^„=(° S/i? ). (A5, 
Here g^(x) , g^(x) - the orthonormal eigenfunctions of the operator d 2 x : 

d 2 x gi = -\ A m gi , dlg B m = -^ n g B m (A6) 

They are supposed to be defined at the segment [— L, L] and to satisfy the boundary condi- 
tions for 771 , r/2 respectively. 

The functions /^(y), fn(y) are the orthonormal wave functions of harmonic oscilators 
satisfying the conditions A + Af^(y) = nf^(y) , B + Bfj/{y) = nf^(y). They may be repre- 
sented as follows 

1 y2 1 u 

= 71-1/4(1 - Z/) 1 ^ 6 2<1 ~" )1/2 ' = ' f» = ( n !)V2 (^ + ) ^ 



/o B = ^(T^ e "^ • *'' = ° • t'Wfit. 6 *)" 1 *' (A8) 

For the perturbation operator one has 



= 0, = 0, = -2e(g**d x g£)(Cd y tf 



1 J 
(A9) 
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Due to relations (A. 5), (A. 9), the expression for the minimal value of a is a = 
min{E A , E B }, where the quantities E A ,E B are reduced in the second order of perturba- 
tion theory in the operator V to the form 



F A_ F A y 45 2 | < gg\d x \g£ > | 2 | < fjMtf > 



m,n 



pB _ f b y 4e 2 | < gj\d x \g B > | 2 | < fMfl > | 2 

m,n ran 00 

Eigenvalues of the unpertubed Hamiltonian H may be written as 

E A n = (l + u)\ A l + 2VT^(n+l) 

E B n = (l-u)\ B l + 2VTT^(n + l) 
Inserting (A. 7), (A. 8) into the matrix element < fn\dy\fo >, one gets 



(A10) 
(All) 



(A12) 



<fn\d y \fo B >= J ffd y f B dy = 



-i n (l - z/)™ /4 



( n |)l/2 2 n/2 7r l/2( 1 _ ^2)1/8(2 + ^1/2 

The relation 



ye 



OO , . m — 

J = / dye 2 y[ay - d y ) e 2 = V27r2 2 a 2 ^ -f__ y;j! i 



permits to transform (A. 13) to 

< fn\dy\f B >= < 







i n V2(l - z/) 3/8 (l + z/) 1/8 (n!!)C ra/2 (s^) (n - 1)/2 

(n!) 1 / 2 ( v / TT^ + ^T^) 172 !^! 1 / 2 



where 



C 



VI + f- 



(A13) 

, m = 1, 3 . . . 

(A14) 

n = 0,2,4,... 
n= 1,3,5,..., 

(A15) 

(A16) 



Vl + ^ + VT^V 

The expression for the matrix element < fn\dy\f£ > may be obtained from (A. 15) by 
changing v — > — 1/. 

With the aid of (A. 12), (A. 15) the Eq.(A.ll) is reduced to 



v- 
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£ 2( 1 _ C 2y/2 
(l-v 2 ) 1 / 2 



E 



n=l,3,5,. 







3<f > 






n+i( 





(A17) 



The sum over n is calculated along the following way 

1 i , ..sn.„ i A i 1 



n=lA... n! ^ + A) 



l C n -\n\y)H n+x -^ dt _ r x y, (Ct) n -\n\\) 2 



/or 



n=l,3,... o 
C 



— f — 

'W (1- 



-dt 



A 



arcsin C 



Y / sin A (y)*/. (A18) 



(Ct)2)3/2- ^A+iy (l_ f 2)3/2 jT^C 2 C X+1 

Introducing the designation (3(v) = [1 — (1 — z/) 1 / 2 /(l + z/) 1 / 2 ]/2 and using (A. 18), one 
can transform (A. 17) to the following expression 



£ 2 (l _ C 2 )3 /2 j t P(-») dt 



1-1,2)1/2 y (i_c 2 * 2 ) 3 / 2 



E \<A\g$>?V A M 



(A19) 



m=0,l,2,. 



where 7 ^) = [a£(1 + - i/) - A B ] /2. 

For one gets in analogous way 

1 - C 2 t 2 )3/ 2 m= ^ 



Z/2)V2 



<^l^bo A >|V-M, (A20) 



where = V / T + ^ [A£(l - i/)/(l + i/) - Aq] /2. 

Since one should use now the explicit expressions for functions g A , g^, below we consider 
the cases of specularly reflecting and duffusively scattering boundaries separately. 



a) Specularly reflecting boundaries 

For the given crystal orientation the boundary conditions for functions g A , g B in the 
case of specularly reflection are 

<?H*=±l = , d x g* lx=±L = (A21) 

The orthonormal eigenf unctions of the operator d 2 satisfying the conditions (A. 21), consti- 
tute the sets 

9 ™ = ~7l C0S ( km+lX + ^ m+1 ) ' 9 ™ = ~j=^ sin(k m x + (p m ) (A22) 
where k m = jr-m , <p m = |(m + 1) , m — 0,1,2,... 
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Substituting (A. 22) into (A. 19), (A. 20) and taking into account relations 

, m + 1 
-k\ , m — 1 



<gt\d x \g*>=0 , <&\d x \g*>=\ C I' Zt ! (^-23) 



one obtains 

E B = VlTT> (A24) 



Eq. (A. 25), which is applicable for small enough e and any admissible value of v, is reduced 
in the limit \v\ C L 2 to the following expression 

. /ir\ 2 / <r 2 (l-C 2 ) 3 / 2 } \ , 

H «L 2 : ^ = v^ + y b + ,- ( \_^ M 

b) Diffusively scattering boundaries 

Under the boundary conditions 

^|x=±l = , ^|,= ± l = (A27) 
the corresponding eigenfunctions of the operator d 2 x may be written as 

9m = 9m = ^ cos(k m+1 x + <p m+1 ) , (A28) 

where k m = mii/2L , (fi m — (m + 1)tt/2 , m = 0, 1, 2, 

From (A. 28) one gets the following expressions for matrix elements 

1 f 

< g m \d x \go >= — J cos(k m+1 x + (f> m+1 )d x cosfax + cp^dx = 

-L 

2h (1 - (zrm + m = 0, 1, 2, . . . (A29) 

2 



7r \m + 2 mj 



and for the eigenvalues of the operator <9 2 : Xq = kf = , = = (m + l) 2 . 

With the aid of (A.29) the relation (A. 19) takes the form 

^vT+^+(i-,)(^) 2 - 
(i - „ 2 )v 2 / (i - cvfn ti v ^ _ ij 2 i AM) 
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Since g A = g^, the expression for E A may be obtained from E B simply by substitution 
v — > — v. It is easy to see that the curves E B (L) and E A (L) cross each other. Designating 
v = \v\sgn (T — 7*.), where is the point of crossing, one obtains for a = min 



a = Vl + i>+ (1 - v) ( 



7T \ 2 



- - 



2L/ 



/ , f _^_ l( (ft)"^((«)^.) (.4.31) 

L 2(l _ ,,2)1/2 7 (1 _ C2i 2 )3 /2 2. ^ \ ) 

Eq.(A.31) can be simplified in the limits of large and small values of parameter L. For 
L 2 > 1 it follows 

2 



: «=^ + (2l) (1-^ (A32) 

where the relation 



~ A; 2 vr 2 

has been taken into account. In the opposite limit L 2 < 1 the integral in (A. 31) is trans- 
formed as follows 



1 



i 



_ / dt - t 5(i>) = - / — - t m = 

L 2 J (1 - CH 2 fl 2 n 2 VT^ (B 4fc2 " X ) / (1 " CH2f ' 2 



8 1 

7T 2 V^^ (B 4 ^ 2 - l) ( X " C ' 2 ) 3/2 

Here 



(A33) 



Then one has 



7T ^ 2 



L 2 «l : a= ^— ) (l-z>) + v / TT^ 



2e 2 ~ A; 2 
~~ ZJn i ,~,\2 2.^1 



(A35) 

The sum presented here is 

y fe 2 1 tt 2 / 8 7 cot(7r 7 ) \ 

fe =i(A; 2 -l) 2 (^ 2 -7 2 ) l-4 7 2 l vr(l-4 7 2 )J l ' J 
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From (A.35), (A.36) it follows 

(A37) 
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FIGURES 

FIG. 1. The orientations of crystalline axes and the magnetic field relative to the film bound- 
aries. 

FIG. 2. The dependence H C 2(T) for a film from superconductor with two-component order 
parameter for v = e = 0.5 and film thickness 2d = 20£o in the case of speculary reflecting surfaces 
(1) and of diffusive electron scattering (2). 
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TABLES 



TABLE I. 



Singlet 


Triplet 




1 


A lu 


CPz 


Mg 


Im^ + ipyf 


A 2u 


cp z lm(p x + ip y f 


Big 


p z lm(p x + ip y ) 3 


B\ u 


c\m(p x + ipyf 


B 2 g 


p z Re(p x + ip y f 


B 2u 


cRe(p x + ip y )' s 


Elg 


PzPx,PzPy 


E\ u 


CPxi CPy 


E 2 g 


2PxPy,pt-Py 


E 2u 


C2p x PyP Z ,Cp Z (pl-pl) 



TABLE 2. 



Representation 


(()q = nn/3 


O = (2n + 1)tt/6 


Alg, Alu 


rf(0) = 


V(o) = 


A 2 g, A 2u 


?7(0) = 


7/(0) = 


B\g, Blu 


t/(0) = 


7/(0) = 


B 2 g, B 2u 


77(0) = 


7/(0) = 


Elg, E lu 


m(o) = 

77 2 (0) = 


m(o) = 

7%(0) = 


E 2 g, E 2u 


?7i(0) = 
fj' 2 (0) = 


?7i(0) = 

772(0) = 
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